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Abstract 

We study the chiral logarithms in AS = 1 kaon decay amplitudes from new flavor physics in 
beyond-standard-model theories. We systematically classify the chiral structures of dimension-5, 
6 and 7 effective QCD operators constructed out of light-quark (up, down and strange) and gluon 
fields. Using the standard chiral perturbation theory, we calculate the leading chiral-logarithms 
associated with these operators. The result is useful for lattice calculations of the QCD matrix 
elements in K — > tttt decay necessary, for example, to understand the physical origin of the direct 
CP violation parameter e'. As a concrete example, we consider the new operators present in 
minimal left-right symmetric models. 



I. INTRODUCTION 



Non-leptonic kaon decay has been a focus for both theoretical and experimental physics 
for over 40 years since the discovery of CP violation by Christenson, Cronin, Fitch and 
Turlay pj in Kl — > 2ir. Since then the origin of CP violation has long been a challenge to 
many theoretical models. The recent data from various experiments have yielded a clear 
non- vanishing direct CP-violation parameter 0, 0] : 

Re(e'/e) = (16.7 ± 2.6) x 10" 4 , (1) 

which ruled out the so-called superweak theory where no direct CP violation appears in 
the decay jij. At present, a full theoretic explanation to the origin of this phenomenon 
is still lacking. In the framework of the standard model (SM), direct CP violation can be 
generated by the non-zero phase in the quark flavor-mixing matrix (CKM matrix), as was 
suggested by Kobayashi and Maskawa (EJ. A precision calculation of the effect, however, is 
extremely hard due to the non-perturbative nature of the strong interactions at low energy. 
Results from several groups utilizing different methods differ widely, with error bars much 
larger than that of the experimental result 0] . The unsatisfying situation of the theoretical 
calculations have attracted much interest in attributing part of the phenomenon to physics 
beyond SM. 

To be able to pin down the contribution to e' from models containing new physics, one 
has to make precision calculations of the strong-interaction physics associated with the non- 
perturbative structure of kaons and pions. Various methods have been used to calculate the 
hadronic matrix elements, such as lattice QCD-inspired models [9|,[l0j, chiral expansion 



together with large- N c [11|], and parametrizations [12|. At present, the lattice field theory 
is the only approach based on first principles, with controllable systematic errors. There 
are difficulties in lattice calculations which are associated with the fact that the final state 



contains more than one particle. By Maiani- Testa theorem [13j, it is impossible to extract 
the physical kaon decay matrix elements by taking the limit r — > oo in the Euclidean space. 
In practise, there are several ways to avoid it: one can either work with an unphysical 



choice of momenta [141 . Il5j . utilize an unphysical set of meson masses |16l . 117] . or derive 
the physical matrix elements by unphysical, but calculable ones. All of these methods need 
chiral perturbation theory (ChPT). 

ChPT assumes an approximate chiral symmetry exists in SU(3)f and describes the low- 
energy QCD physics under a chiral breaking scale A x ~ m p by the pseudo-Goldstone par- 
ticles, namely pions, kaons and eta. Then the low-energy physics can be perturbatively 
expanded in powers of the particles' external momenta and masses. It further assumes that, 
the Wilson coefficients of the QCD operators expanded in terms of meson operators are inde- 
pendent of the external states. Therefore the amplitudes of a large number of reactions can 
be determined by a relatively small set of coefficients, which gives us the predicting power. 
In the case of kaon decay, ChPT is used to connect the desired matrix element (irn\0\K) 
with some unphysical quantities, such as (ir\0\K) and (0\O\K). The results in ChPT are 
needed before doing relevant lattice calculations. Here we will neglect some subtleties in 
the ChPT (such as quadratic divergence cancelations, zero pion mass corrections, etc) and 
focus on possible operator structures as well as their chiral logarithm corrections for the 
kaon decay process. It is the goal of this paper to examine the chiral structures of possible 
QCD operators responsible for As = 1, Ac? = —1 decay in generic beyong-SM theories and 
to calculate the large chiral logarithms associated with them. Previous calculations have 
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20, 21, 22, 23 . Our work extends these 



been made for operators present in the SM [18|, \li 
studies to all possible operators in new physics models. 

The paper is organized as follows: we start from the operator basis in SM for the kaon 
decay, as well as possible new operators coming from physics beyond SM. A chiral pertur- 
bation theory calculation will be presented in the following section, with all corresponding 
operators and their one-loop corrections of the matrix elements. We end this section by 
applying our result to a specific example. Concluding remarks and outlook are presented in 
the last section. 



II. EFFECTIVE OPERATORS FROM NEW FLAVOR PHYSICS 

In this section, we consider effective QCD operators contributing to CP- violating K — > nit 
decay in a generic weak-interaction theory. There is an extensive literature on this topic in 



the context of SM [2J, [26| . Our focus is on new operators arising from novel CP- violating 
mechanisms beyond SM. We classify the effective operators in terms of their flavor symmetry 
properties under chiral group SU (3) l x SU (3)r when up, down, and strange quarks are taken 
as light. 

The direct CP violation parameter e' for K — > tttt decay is defined as [2 

s/2 Re A) \ReA ReA 2 ) ' 1 ; 

where 5i is the strong-interaction Tin scattering phase shifts, and Aj is the weak kaon decay 
amplitudes: 

A ie iSl = (7171(1 = 0,2)\(-iH w )\K°) , (3) 

where Tiw is the effective weak-interaction hamiltonian which depends on the underlying 
theory of kaon decay. The small ratio uj = Re^/ReAo ~ 1/22 reflects the well-known 
AI = 1/2 rule. Accurate calculations of e' depend on reliable evaluations of the effective 
QCD operators present in 7i w . Our goal in this paper is to classify these QCD operators 
and study their chiral behavior. 



A. Standard Model Operators 

The standard procedure for calculating e' utilizes an effective field theory approach. The 
physics at high energy (or short distance) can be calculated perturbatively and is included 
in Wilson coefficients. The physics at low-energy scales is included in the effective QCD 
operators composed of light flavor quark fields (u, d, s) and gluon fields. Large QCD radiative 
corrections or large logarithms are resumed by solving renormalization group equations. 
The effective operators responsible for the neutral kaon decay have the flavor quantum 
numbers As = 1, Ad = — 1. In the SM, it is well-known that 7i eS consists of following 10 
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operators |25|, [26J : 



Qi 
Q2 

Q3,5 



(SiUj) V -A(Ujdi)v-A , 
(SiUi)v-A(Ujdj)v-A > 

{Sidi)v-A^(c[j(li)v^A , 



q 




(s i d j ) v _ A ^2(q j q i ) VTA , 




V±A , 





V±A , 



(4) 



where (qq')v±A = qL,Rl^qL,R with qi,R representing the left-(right-) handed quark fields. 
The summation in q is over the light-quark flavors; it, d, s; i and j are color indices; and e q is 
the algebraic charge factor for flavor q. The Qi,2 come from the single M^-boson exchange 
tree diagram, and Q3-6, Q7-10 are derived from one-loop gluon and electro- weak penguin di- 
agrams, respectively. In the SM, (lmA Q /R,eA ) is dominated by the QCD penguin operators, 
whereas (Im^/Re^) receives contribution from the electro- weak penguin operators only, 
because the gluon interaction is flavor-singlet and cannot contribute in A J = 3/2 channel. 

A chiral structure analysis of the above ten operators will be useful if we wish to use 
lattice QCD to calculate the relevant matrix elements in kaon decays. Using (m, n) to 
denote a representation of group SUl(3)x SUr(3), where m and n are the dimensions of 
SU(3) representations, it is then easy to see that Qi,2 and (^9,10 belong to (8, 1) and (27, 1), 
<?3~6 to (8,1), and Q-j^ to (8,8) [26(. ChPT calculations have been made to uncover the 
large logarithms associated with these operators, which in turn help to establish relations of 
different matrix elements useful for lattice QCD calculations. For the reader's convenience, 
we have collected the standard chiral results in the Appendix. 

The question is what is the general chiral structure of all possible weak operators that 
might emerge in theories beyond SM? The rest of this section is devoted to addressing this 
question. 

B. Dimension-5 and -6 Operators 

Let us systematically consider the possible operators and their chiral structures in a 
general low-energy description of kaon decay, independent of the underlying short- distance 
flavor physics that can be taken into account by Wilson coefficients. The lowest dimensional 
operator is a dimension-5 chromo-magnetic operator, 



This operator does appear in the standard model through penguin diagram as shown in 
left-panel Fig. 1, although it is proportional to the strange or down quark masses, which 
is chirally suppressed. It also appears naturally in the left-right symmetric model (LRSM) 
with left-right handed gauge-boson mixing, proportional to charm or top quark masses [27[ , 
as shown on the right panel in Fig.l. Under chiral symmetry, this operator transforms as 




(5) 
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(3,3) + (3,3). The chiral logarithms appearing with the matrix elements of this operator 
has been studied before in the literature, and is collected in the Appendix. 




FIG. 1: Dimension-5 effective operators generated from the weak-interaction vertex corrections in 
SM (left) and in LRSM (right). The crosses on fermion lines represent mass insertion, needed to 
flip the chirality of the quarks. 



Next, consider dimension-6 four-quark operators. We define the following flavor tensor 

&H = (qr qj )(q k T' qi ) J (6) 

where the flavor indices k, I go through 1, 2, and 3, or up, down and strange quarks. T 
and r' are possible Dirac matrix structures. In addition, there are two independent color 
structures (1)(1) and (t a )(t a ) which are not essential for the following discussion. Assuming 
that all fields are projected to their helicity states, the possible helicities are as follows: 

• All four quark fields have the same chiral projection 

• Both q 1 and q k (also j and /) have the opposite chiral projection 

• Both q 1 and q k (also j and /) have the same chiral projection 

These are only possibilities because the operators must be Lorentz scalars, and the numbers 
of left and right-handed fields apart from the first case must be exactly 2, respectively. In 
the first and second case, the operators are the ones appearing in the SM weak interactions, 
as shown in Eq. (J3J), and their parity partners. They correspond to chiral structures (8,1), 
(1,8), (27, 1), and (1,27) from the first case, and (8,1), (1,8) and (8,8) from the second 
case. The new (1,8) and (1,27) structures will appear in, for example, LR symmetric models 
where the right-handed gauge boson plays the same role as left-handed one in the SM. Since 
the strong interactions conserve parity, the new operators in LRSM have the same matrix 
elements as Q^s in SM up to a parity sign. The corresponding Feynman diagrams in both 
SM and LRSM for the first and second cases are shown in Fig. [2J 

In the last case, there are new chiral structures arising from operators of type 

(s L Td R )(q L Tq R ), (s R Td L )(q R Tq L ) . (7) 

where two T's must be the same. The new chiral structures are (6, 6), (6, 3), (3, 6), and their 
parity conjugates. However, (6, 3) and (3, 6) involve symmetrization of two flavor indices 
and, at the same time, anti-symmetrization of the other two. It is easy to check that the 
result vanishes, and we are left with just (6, 6) and (6, 6). 
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FIG. 2: Feynman diagrams generating dimension-6 quark operators in SM and LRSM. 



Let us consider the following tensor with up (and hence lower) indices symmetrized 

©&) = ItfrPqRkMrqm) + \{&<lmWi?<lm) , (8) 

where the upper indices represents the left-handed fields and the lower indices the right- 
handed. Without loss of generality, we take % = 3. If j = 3, and k or / is 3 and the other 
indices must be a 2, one get an isospin-1/2 operator 

e VM = 33 = s L Td R s L Ts R . (9) 

If we define a tensor T^j which multiplies the quark operator Ofj to generate the above 
operator, T^Ofj, we have 

if = if = 1/2 , (10) 

and other components zero. 

On the other hand, if j, k, I take l's and 2's, one can subtract the trace with respect to 
j and k, and j and /, and one obtain an isospin-3/2 operator 

0(6,6) = R 12 , R 21 _ ft 22 
U 3/2 — U (31) + U (31) U (32) 

= s L Tu R u L Td R + s L Td R u L Tu R - s L Td R d L Td R , (11) 
with corresponding non-zero tensor components 

t~>31 t->31 t->13 ml3 t->32 t->23 I/O (~\0\ 

1 12 — 1 21 — 1 12 — 1 21 — 1 22 — 1 22 ~ i / z ■ ( iz J 

Another isospin-1/2 operator will can be obtained its trace part, 

R (6,6) ft l2 , 021 , oo22 

U l/2,5 — U (31) + U (31) + ZU (32) 

= s L Tu R u L Td R + s L Td R u L Tu R + 2s L Td R d L Td R , (13) 

and the corresponding tensor components are, 

Til = Til = Til = T 2 \ 3 = 1/2; if = if = 1 . (14) 

An example of these new operators in LRSM through flavor-changing neutral and charged 
currents is shown in Fig. [3l The relevant QCD four-quark operator will be: 

As=1 = fad*) 

1 

= \ K? S + - e<«] (15) 

where 6 (3 ' 3) = (s L u R )(u L d R ) - (s L d R )(u L u R ). 
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FIG. 3: Feynman diagrams generating scalar quark interactions through neutral-current Higgs 
exchanges. 



C. Dimension-7 Operators 

Dimension-7 operators come in two types. The first is the chromomagnetic operators with 
an insertion of two additional derivatives, which does not change the original chiral structure. 
The second type is an insertion of one derivative into four-quark operators discussed above. 
Since the covariant derivative has one Lorentz index, it must be contracted with another 
Lorentz index appearing on a Dirac matrix. An example of this type of operators is 

o L P R g = hd^ndj,Rj2^ L ^ D ^ L > 

g 

£>p R ew = hL{^ V )di,R^2e q q jtL { lfl D v )q j!L , (16) 

g 

which comes from the gluon and electromagnetic penguin diagram in LRSM as shown in 
Fig. [TJ These operators contain either 3 left-handed fields and 1 right-handed one, or 3 
right-handed fields and 1 left handed one. They have novel chiral structures (15, 3), (15, 3), 
(6, 3), (6, 3), and parity partners. Let us classify them all in details. 

1. (15,3) 

Let us use 0^- to represent an operator ^rg ife g^r'g ra , where I is flavor index of the 
right-handed field, and V and V are not just Dirac matrices. We first construct 0fj which 
forms (15, 3) after symmetrizing the up two indices and subtracting the traces, 

2eg = eg + e« - 1% K- + e S] - \% K + ©a > ( 17 ) 

where a sums over 1, 2, and 3. Clearly, i, or equivalently j, has to be an s. One isospin-1/2 
operator that one can immediately identity is when j is 3, k is 2, and I is 3, namely 

e;f 2 3) = el = e*j = s L vd L s L v's R , (is) 

where the only non-zero tensor component is 

Til = 1 . (19) 



7 



Other independent operators can be obtained by considering j, k, and / as up and down 
quarks. Others, such as ©33, can be related to these through traceless conditions. 

One can get an isospin-3/2 operator by symmetrizing k and I while taking away SU(2) 
traces between j and k, and j and /. Thus we have the following combination, 

©S + ©8 - l*j (®t + ©k) - 1$ (©3« + ©3a) , (20) 

where a sums over 1 and 2 only, and j, k and I can take value in 1 or 2. There is only one 
independent operator 

e 3 ™ = 2 (e- + ell - ©1 

= s L Yu L u L Y'd R + s L Td L u L T'u R + u L Tu L s L T'd R 

+u L Td L s L T'u R - s L Td L d L T'd R - d L Td L s L T'd R . (21) 

The corresponding tensor components are 

T~~>31 rpl3 rp31 rpl3 rp32 rp2'3 1 ( r ) r )\ 

1 12 — 1 12 — 1 2l — 1 21 — ~ 1 22 — ~ 1 22 ~ 1 • \ ZZ ) 

The trace part of the above operator produces an I = 1/2 operator, ©^ + ©3". We can 
subtract from the result with another isospin-1/2 operator, ©I 3 ,, to cancel the unwanted 
trace part — |5f [©"'•] ~ \^ [©«'] m Eq. ffTTl) . The resulting / = 1/2 operator is 

©SS s 2 (©3? + ©" + 2©l) " 3©| 

= s L Yu L u L Y'd R + u L Tu L s L T'd R + s L Yd L u L Y'u R 

+u L Td L s L T'u R + 2s L Td L d L T'd R + 2d L Td L s L T'd R - 3s L Ts L s L T'd R , (23) 

with the following tensor components, 

T^31 _ 7U3 _ o-i31 _ 7-1I3 _ 1t^32 _ 1t^23 _ 1^33 _ i / 9z i\ 
J 21 ~~ 1 21 ~ 1 12 ~ 1 12 ~ 2 22 _ 2 ~~ 3 2 _ 



Finally, one can antisymmetrize k and / to generator another isospin 1/2 operator. Again 

32 
33 



we add ©33 to cancel the unwanted trace part: 



©iS - 2 (©31 - ©"J + ©33 

= s L Td L u L T'u R + u L Td L s L T'u R - s L Tu L u L T'd R 

-u L Tu L s L T'd R + s L Ts L s L T'd R , (25) 

with the following tensor components 

Til = T 2 \ 3 = —T^l = -Til = Til = 1 . (26) 

Note that the operators in (3, 15) can be obtained from the above through parity transfor- 
mation. 
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2. (3, 15) 

Define 0fj = qjXq R k(? R T'q R i, an( ^ construct the general (3, 15) operators, 

2eg = eg + eg - U k 3 [eg + eg - -5} [e« fc + eg , (27) 

where the a trace is over 1, 2, and 3. Either index % or j can be identified as the strange 
quark field. In either case, one can construct isospin-3/2 operators by subtracting SU(2) 
trace. With the left-handed strange quark, we have 

- ©a? + eg - el = e\l + eg - e| 

= s L Tu R u R T'd R + s L Td R u R T'u R - s L Td R d R T'd R , (28) 
with tensor components, 

Til = Til = -Til = 1 • (29) 
With the right-handed strange quark, 

B (3,15) A12 , A21 _ A22 _ Q 12 , Q 21 _ Q 22 

W 3/2,i? — w 13 ^ w 13 w 23 ~~ w 13 ^ w 13 w 23 

= u L Tu R s R T'd R + u L Td R s R T'u R - d L Td R s R T'd R , (30) 
with tensor components, 

Tll = TH = -Tll = l. (31) 
There are also two isospin-1/2 operators. The first one with left-handed strange quark, 

- 4 (eg + el\ + 2 eg) 

= s L Tu R u R T'd R + s L Td R u R T'u R + 2s L Td R d R T'd R 

-3s L Ts R s R T'd R - 3s L Td R s R T's R (32) 

with tensor components 

1 1 1 

rp31 rp3l _T~'32 _T~'^^ — 1 

1 12 ~ 1 21 ~ 2 22 ~ 3 32 ~~ 3 _ 

And the second one has the right-handed strange quark, 

©53 = ^ + ©?3 + 2ei = eg + eg + 2 e* 

= u L Tu R s R T'd R + u L Td R s R T'u R + 2d L Td R s R T'd R (34) 

with tensor components 

^i2 3 = T 2 \ 2 = ^11 = 1 • (35) 
The efj is not independent by the same reason as for (15, 3). 
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3. (6,3) 

Define 0f- = q l L Tqi,k(? I T'qm, an d construct the (6,3) operator, 

U ij — e «jmV — r^m \ Z V a/3 + e U Q/3 j . {60) 

where k, m are symmetric, and a, /? run over 1 to 3. To get the As = — Ad = 1 operators, 
none of the m, k and I can be a 3: when 3 is on the e, it prevents both % and j from being 
a strange quark, and when 3 is a lower index, both i and j must be 3 which is impossible 
because of the antisymmetry. In fact, the only possible combination for m, k and I is 2, 2, 
1. 

To get an isospin 3/2 operator, one must symmetrize m, k and I, yielding 

ef f f = e 221 + e 212 + e 122 = el\ + + e 22 

= s L Yd L u L Y'u R + s L Yu L u L Y'd R - u L Yd L s L Y'u R 

-u L Yu L s L Y'd R - s L Yd L d L Y'd R + d L Yd L s L Y'd R . (37) 

The corresponding tensor components are 

rpZl T^13 T^31 rpl3 T^32 rp23 1 

1 12 — L \2 — L 2\ — 1 2l — J 22 — J 22 — 1 • l oc 7 

There is also an isospin-1/2 operator by anti-symmetrizing k and I 

ef /2 3) = 29 221 - e 212 - e 122 = 2 (e 2 } - 

= 2s L Yd L u L Y'u R - 2u L Yd L s L Y'u R - s L Yu L u L Y'd R 

+u L Yu L s L Y'd R - d L Yd L s L Y'd R + s L Yd L d L Y'd R , (39) 

The tensor components are 

l -TH = = -Til = Til = Til = -Til = 1 • (40) 

4- (3,6) 

Define 0fj = q l L Yq Rk q R Y'q R i, and construct the (3,6) operator, 

eg = e klm Q i{jm) = -e klm (e a/ 3 m Qff + e a/3j Qf£j , (41) 

where j, m are symmetric, and a, (3 run from 1 to 3. The only choices which will generate 
the As = —Ad = 1 operators is i, j and m takes 3, 3, 1. If i is 3, and j and m take 3 and 1, 
we have the isospin 1/2 operator 

®1/2,L = ®313 + ©331 = 203! 

= s L Yu R u R Y'd R - s L Yd R u R Y'u R + s L Yd R s R Y's R - s L Ys R s R Y'd R , (42) 
with the following tensor components, 

rp31 rp31 T"i33 T">33 1 /A'i} 

1 12 — 1 2l — 1 23 ~ 1 32 — L ■ \^°) 
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On the other hand, if % takes 1, one gets another isospin 1/2 operator 

ft(3.6) _ A _ A12 
U l/2,i? = U 133 - U 13 

= u L Tu R s R T'd R - u L Td R s R T'u R , (44) 

with the following tensor components 

TH = -Tg = 1 . (45) 

One could consider operators with dimension 8 and higher. However, generally they are 
suppressed by 1/A 2 relative to those we have considered, where A is some weak interaction 
scale. We summarize the above result in the following table: 



TABLE I: Chiral representations appearing in dimension-5, 6 and 7 operators. 





dimension-5 


dimension-6 


dimension- 7 


(L,R) 


(3, 3), (3, 3) 


(8,1), (1,8) 
(27,1), (1,27) 
(8,8) 
(6, 6), (6, 6) 
(3, 3), (3, 3) 


(15, 3), (3, 15), (15,3), (3, 15) 
(6, 3), (3, 6), (6, 3), (3, 6) 
(3, 3), (3, 3) 



We emphasize that these operators are completely general, independent of the underlying 
mechanisms (supersymmetry, large-extra dimension, or little Higgs, etc) for flavor and CP- 
violations in beyond SM theories. 



III. CHIRAL EXPANSION AT LEADING ORDER 



Chiral perturbation theory (ChPT) for kaon decay is useful for two reasons: First, it al- 
lows one to connect the physical matrix elements (K\0\im) to some unphysical, but easier- 
to-calculate matrix elements on lattice. Second, it yields dependence of the matrix elements 
on meson masse parameters. Since lattice calculations are usually done at larger and un- 
physical meson masses because of limited computational resources, this dependence can be 
used to extrapolate the calculated matrix elements to physical ones. In this section, we will 
build a set of effective operators in ChPT up to the lowest order, and use an example to 
illustrate how to connect the unphysical processes to the physical process K — > tttc we are 
interested in. 

In lattice calculations the quenched approximation to QCD has usually been applied in 
the past, where valance quark fields are "quenched" by corresponding ghost quark fields with 
the same masses and quantum numbers but opposite statistics. The ChPT can be adapted 
with the quenched QCD by introducing the "super-?]"' field into the effective lagrangian (28| . 
In this paper we will work with the full dynamical QCD only. 



A. ChPT and SM Operators 

The standard ChPT starts with the nonlinear Goldstone meson field S by: 

f 



E = exp ( ^ ) , (46) 
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where d> is the Goldstone meson matrix 



' V2 + Ve 



\ 



7T 

K~ 



7T 

#0 



'/ 

V6 



K+ \ 
K° 

fn j 



(47) 



and / ~ 135 MeV is the bare pion decay constant. We separate the effective ChPT La- 
grangian into two parts: 



C 



ChPT 



£a ~\~ (Li 



(48) 



where C s corresponds the QCD strong interaction which preserves the flavor symmetry; the 
C w is an effective Lagrangian for non-leptonic weak interaction, and is responsible for the 
As = 1 processes. The lowest-order terms for the strong interaction part is: 



f 2 

£(2) = ^_Tr(^S^S) + vTr [ME + (ME) 1 '] 



(49) 



where M = diag(m ?t , ma, m s ) is the quark mass matrix; and v ~ —\{uu) is proportional to 
the quark chiral condensate at chiral limit. We demand the fields transform under SU(3)l x 
SU(3)r as: 

S -> LY>R\ M -> REL\ (50) 

to keep the Lagrangian invariant under an SU(3)l x SU(3)r transformation. Higher order 
terms in the effective Lagrangian contain higher derivatives, and can be written in system- 
atic derivative and mass expansion. For our purpose here, however, only the leading large 
logarithms are calculated, and the higher order terms are irrelevant. 

At one loop, the physical masses and wave- function renormalizations are given by [2J 



m 



7T.0 



rn 



K 



rn 



K,0 



1 + L(m n ) - -L(m v ) + ... 



1 + -L(m v ) + 



4 2 
1 + -L(m^) + -L(m K ) + 



Z K 

u 

fx 



1 + -L(m fl 



^L(m K ) + ^L(m v ) + ... 



f [1 - 2L(m T ) - L(m K ) 



f 



-L{m n ) - -L(m K ) - -L(m v ) + 



4 v n, 2 v ,w 4 
for the pion and kaon fields, respectively. L(m) is the chiral logarithm defined as: 



(51) 
(52) 

(53) 

(54) 
(55) 

(56) 



L(m) 



i i 
m m 



(57) 



with fi x the cutoff scale. The dots represent non-logarithm contributions from 0(p A ) and 
higher order Lagrangian terms. In this paper we focus only on the large chiral logarithmic 
corrections and will not, for simplicity, include the dots explicitly in the results. 
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In the standard electroweak theory, there are 7 independent four-quarks operators which 

n 



can be classified into (8, 1), (27, 1), and (8, 8). Define 6 = T^q^^^q^q^i^qu-, we can obtain 
four-independent quark operators with the following tensor components 



(27, l) 3/2 


m31 
J 21 


m31 
— 1 12 


rpl3 

— 1 21 ~ 


rn\3 7^32 
1 12 ~ 1 22 ~ 


7^23 
J 22 " 


1 

= 2 ' 




(27, l) 1/2 


7^31 
1 2l 


7^31 
— 1 12 


rpl3 

— 1 21 ~ 


7^13 1 7^32 
J 12 — 2 2 ~~ 


7"i23 
J 22 — 


1 7^33 7^33 
L i 1 23 ~~ J 32 


3 

~~ ~2 


^8, l)l/2,5 


rp3l 

1 21 


th13 
_ J 12 


rp31 

— 1 12 — 


7^13 1 7^32 
J 21 — 2 2 ~~ 


7^23 
J 22 — 


7->33 7133 
J 23 — J 32 ~~ 


1 • 


8, l)l/2,A 


rp3l 

1 2l 


7^13 
— J 12 


rp31 

— 1 12 


7^13 1 

- J 2i - 2 • 









(58) 
(59) 

(60) 



On the other hand, defining a (8,8) operator = T^q^^^q^q^^q^i, we have three quark 
operators with following tensor components, 

(61) 
(62) 
(63) 

One can similarly defined other (8,8) operators with a different color indices contractions. 
For the sake of convenience, we have broken the operators into representations of definite 
isospins. This has the advantage of easily building up reducible operators from linear com- 
binations of these simple ones. For example, the SM electromagnetic penguin operators 



3)3/2 


. 7^31 
• 1 21 


7->31 
— 1 12 ~ 


rp32 -1 
1 22 — 1 


j 


1/2,5" 


. rp31 

■ 1 2l 


7^31 
_ J 12 _ 


T£/2 = - 


-T 2 f/3 = 1 


1/2,A 


. 7^31 
• J 21 


7^31 
_ 1 12 


rp33 

— 1 2S — 


1 . 



Q- 



e 



1(8,8) 



(64) 



4- V 
2 [ 3/2 ^ W 1/2,A 

and Qs is similar to Q7 but with different color indices contraction. 

In ChPT, one can match the above QCD operators to the hadronic operators made of 



Goldstone boson fields [23|, |30|, [31 

»(8>l) 



ei 8,1) = Tr [A<9 At E<9 /i E t ] 



e 



(8,1) 







(27,1) 
Al 



0(3,3) 



— Tr [ASM + A(SM)t] , 
Tr[AE f ] , 



(65) 



where A = <5j, 3 <5j,2 and T's are tensor structures defined above. The expansions go like 







(8,1) 



a 



(8.1)A(B,1) 



©r ; +«2' 1) 02 8,1) + 



I? 







(27,1) 



Al 



a* 2 ™ eg™ 



Al 



»0 



Al 



+ 



(66) 



where a^ L ' R ^s are "Wilson coefficients" which are universal in different processes, and dots 
represent higher dimensional operators. The subscript i on the (8, 1) operator indicates 
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different quark operators in the same chiral representation, including ones with two right- 
handed fields coupled to the singlet. 



The one loop results of these operators in various processes can be found in [18|, ll9j, [2 



21 . 22 . 23[. Due to the different definitions of the operators and the nonlinear meson fields, 
there might be sign differences among these results. Overall speaking, the (8,1) and (27,1) 
operators dominate in the CP-conserve process. The (8,8) operators, corresponding to the 

1,S3- 



Q 7t 8 operators, play a significant role in CP- violation processes 



The lowest-order 

mass-dependent term OJf'^ will vanish in physical process K — > hit to all orders. This 
property was pointed out by jUJ first and has been well-studied by jH, 
back to this issue later. 



34 . We will come 



B. Chiral Matching of New Operators 

Now we can proceed in constructing new hadronic operators for new interactions arising 
from physics beyond SM. We label operators by the irreducible representatives and their 
isospin quantum numbers. Similar to the case in SM, we define the effective operators at 
their lowest order as: 



Dimension- 6: 
Dimension-7: 







(6,6) 
AI 



.i ■ 



e 



(15,3) 



Al 



AI j^St)*^)!. 



2>(3> 15 ) 







AI 

a (3,6) 







Al 



p-,(15,3)-|j; 

[Tf; 15) ]^(^st)^(st9 M s)^. 

p-i(3,6) yij 



AI IkH 

^Et)*(£t^£V. 



We can also construct operators with one insertion of quark masses, 

Xf = (SM) ± (SM) f , Xf = (MS) ± (MS) f . 



They transform under SU(3)l x SU(3)r as: 

X£ -> LX±L^; Xf 



(67) 

(68) 
(69) 
(70) 
(71) 

(72) 

(73) 



With the insertion of X± we can build two additional sets of the dimension-7 operators at 
the lowest order: 



U AI,X ± 



[T 



for (L,R) belong to (15,3) or (6,3), and 



(L,Rhij 



Al 



J A/ 



(X: 



U AI,X ± 



±) 3 ' 



(74) 



(75) 



for (L,R) belong to (3,15) or (3,6). Therefore, the dimension-6 and -7 QCD operators 
should be matched to hadron operators as follows, 



^D6 
^D7 



a (L,R)Q(L,R) ? 

a (L,R)Q(L,R) + a ^fQff) + a%*>&Vf> , 



(76) 
(77) 
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where higher-order terms have been omitted. 

In SM we need operators with X + only since all QCD operators obey the CPS symmetry, 
the CP transformation followed by an exchange of s and d quarks [3l| . However, the 4-quark 
operators derived from new physics do not necessarily have this symmetry, and hence we 
can have an additional set of operators in the effective theory. 

Just as the (8, 8) operators in SM, the (6, 6) dimension-6 operators will contribute at 
O(p ) order in ChPT. This set of operators can be derived from the Higgs (or some new 
heavy bosons) exchange and will contribute to the CP violation phase in the same manner 
as Q 7 , Qg in SM. We will consider an example of applying our result later. 



TABLE II: Tree level contributions from dimension-6 operators 



(L,R)ai 


K° — > Vacuum 


K+ -> 7T + 


K° 7T ^ 


bo 


CO 


d 


(6, 6)3/2 





-4 


-8 


(6, 6)i/ 2 ,s 


-6 


-10 


16 


(6, 6)i/ 2i A 


-2 


-2 






When calculating CP conserving matrix elements, the new operators are usually negligible 
compared to the SM weak-interaction operators defined in Eqs. fl58|) - fl65|) . The new operators 
are ma inly responsible for the CP-violating phase, and are worth investigating as the case 



111 



321 ] . We observe that some of these new operators, notably (6,6), (15,3), (6,3) and 
(3,15) have contributions in AI = 3/2 channel, in addition to (27,1) and (8,8) operators 
in SM. Furthermore, as we shall see, they all receive large chiral logarithmic corrections 
in one-loop ChPT. Therefore, the operators from new flavor theories beyond SM can help 
explain the AI = 1/2 selection rule and the direct CP violation parameter e'. 



C. Results at Tree Level 

In this subsection, we consider tree-level relations among the matrix elements of the QCD 
operators in different states. These relations reflect chiral symmetry and can also be derived 
using old-fashioned current algebra. 

There are three processes that we are mainly interested in: K° — > Vacuum, K + — > tt + 
and K° — > 7r°7r°. For K° — > 7r + 7r~, one can obtained the matrix elements through angular 
momentum relation, as shown in Appendix. At tree level, the dimension-7 momentum 
operators will not contribute to K° — > Vacuum process. For the mass-dependent operators, 
the result will be either proportional to (m s — md) ~ m2 Ko — m "i or ( m « + m d) ~ m |-o> by 
the lowest order expansion of X±. However, beyond tree level, the result will no longer be 
proportional to (m s — m d ) or (m s + m d ). 
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TABLE III: Tree level contributions from dimension-7 operators 



(L,R)ai 




> Vacuum 




K + 7T+ 




K° 


^7r°7T 




6 0,+ 


K- 


Co 


c o,+ 


c 'o,~ 


d' 


d o,+ 


d'o,- 


(15, 3)3/2 








-8 





2 


8 


-2 


-2 


(15, 3)!/ 2 ,5 


9/2 


3/2 


-8 


-3/2 


2 


8 


-2 


-2 


(15, 3)l/9 A 


-1/2 


1/2 





-1/2 





8 


-2 


-2 


(15, 3)i/ 2 


1/2 


1/2 





-1/2 














(6, 3)3/2 

V 7 / Of £ 

















8 


-2 


-2 


(6,3)l/2 


3/2 


3/2 


12 


-3/2 


-3 


16 


-4 


-4 


(3, 15) 3 / 2 ,L 








4 





-1 











(3, 15) 3 / 2) fl 








4 





-1 


8 


2 


-2 


(3,15)i/ 2 ,L 


-3/2 


-3/2 


4 


3/2 


-1 











(3,15)i/ 2ij r 


-3/2 


3/2 


4 


-3/2 


-1 


-16 


2 


-2 


(3,6) 1/2 ,l 


3/2 


-1/2 


4 


1/2 


-1 











(3, 6)i/ 2 ,fl 


-1/2 


1/2 


-4 


-1/2 


1 





2 


-2 



For dimension-6 operators in (6,6), we have the tree level results: 

ib 



(0\e D6 \K°) Tvcc = -ja m , (78) 
(n + \e m \K + ) Tlee = ^a m , (79) 



(7r 7r°|e D6 |^°) Tree = l -^a m , (80) 



with bo, Co and do coefficients listed in Table [TTJ which are different from different isospin 
projections. The corresponding results for 7r + 7r~ final state can be obtained from relations 
in Appendix. The non-perturbative coefficient is the same for different operators and 
final states. 

Similarly, the tree-level matrix elements for dimension-7 operators are: 

4277 

(O|0D 7 |A'°) Trcc = — [b 0)+ (m s - m d )a D7!+ + b' Q _{m s + m d )a D7 -] 

= j [b ,+ ( m K,o- m lfi) a D7,+ + b' _m 2 Kfi a D7 _] , (81) 

9 

777 

(ir + \<~) D7 \K + ) Tmc = — j^-[c' a D7 + c' 0+ a D7i+ + cd_a D7 _] , (82) 

9 

777 

(7r 7r°|e I)7 |J^ )T re e = -^{d'o^D 7 + d' Q+ a D7)+ + (_« D7r ] , (83) 

with coefficients listed in Table [mi Here the non-perturbative coefficients a's are different for 
different chiral representations. From the above equations, it is clear that one can obtained 
the two-pion matrix elements from the vacuum and one-pion ones, if only one of the X + 
and X_ types of operators is present, such as in the SM case. 



16 



IV. CHIRAL LOGARITHMS AT ONE-LOOP 



ChPT calculations of the kaon-decay matrix elements up to higher chiral orders are 
needed for understanding the size of chiral corrections and for extrapolating matrix elements 
from unphysical quark masses to physical ones. In lattice calculations, unphysically large 
quark masses are usually used to make calculations feasible. Then one needs to extrapolate 
the matrix elements to the physical region. In this section, we calculate the large chiral 
logarithms of the dimension-6 and 7 operators for the process K — > 0, K — > n and K — > tttt 
hoping to get the leading corrections as the function of quark masses. 

In our calculations, we have made the simplifying assumption m u = m&. For (0\O\K°) 
matrix element, we have kept all the Goldstone boson masses independent. For the matrix 
element (tt + \0\K + ), we utilize a common mass tum for all the mesons to conserve momen- 
tum. In the calculation of (7r ir \O\K ) matrix elements the pion masses are neglected. Since 
m^/m? K pa 10" 1 with physical pion and kaon masses, this is a reasonable approximation for 
the physical processes. 

A. K° — ► Vacuum 

The diagram we need to consider is shown in Fig. H] below, where and henceforth the 
square dot represents an effective weak interaction operator while the round dots repre- 
sent strong interaction insertions. We have not shown the wave function renormalization 
diagrams, but they have to be included in the final result. 

— ' «^ 

/ \ 
i \ 

\ I 
\ ' 
\ / 

FIG. 4: Feynman diagram for K° — > Vacuum at one loop. 

For dimension-6 operators, the results for (OlOI-fT ) up to one- loop can be written as: 

(0\e D e\K°) = ^[b + b v L(m 2 v ) + b K L(m 2 K ) + KL{mD] , (84) 

with coefficients listed in Table IIVI Note that we have three different chiral logarithms 
corresponding to eta, kaon, and pion, respectively. The isospin-3/2 operator does contribute 
for the obvious reason. 

For dimension-7 operators, the results are more complicated, 

(O|0 D7 |K°) = (0|e D7 |^°>Trcc + y{ [b' v mlL(ml)+b' K mlL(ml) + b' T mlL(ml)]-a D7 

rt,K m2 K 

± 

+ (K tK m 2 K + b' n ^ml)L(ml)]- a D7>± } , (85) 
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TABLE IV: One loop contributions from dimension-6 operators 



(L,R) AI 




— ► Vacuum 






b v 


b K 


K 


cm 




(6, 6)3/2 











112/3 


80/9 


(S,6)i /2| s 


1/2 


33 


57/2 


196/3 


-1024/9 


(6,6)1/2,;! 


25/6 


15 


3/2 


28/3 


-4 



where each chiral logarithms now have different meson mass factors. The coefficients are 
listed in Tables ED and ED 



B. K~ 



7T 



For the K + — > 7r + matrix elements, we utilize a common mass for mesons m \ 



M 



mz 



m 



K 



in the calculation. The Feynman diagrams are shown in Fig. [5j The matrix elements up to 
the leading chiral logarithms are: 



(n + \e m \K + ) 



P 



[co + c M L(m 2 M )] 



(n + \Q m \K + ) = (7r + |e D7 |K + ) TVee 



m 



f 



A ^L(m 2 M ) 



d M a D7 + c' M+ a D7:+ + d M a D 7,- 



(86) 



(87) 



for dimension-6 and 7 operators, respectively. The coefficients are listed in Tables IIVI and 

Em 



FIG. 5: Feynman diagrams for K + — > 7r + at one-loop. 



C. K° vr 7r° 



For K° — > 7r°7T , we neglect the pion mass in the calculation. The diagrams we need to 
consider for are shown in Fig. [6J The results are: 



<7rV|e D6 |;o 



[d + d K L(m 2 K )] 



(7r°7r°|e D7 |K ) = (7r°7r°|e D7 |K%ee 



m 



f 



d' M a>D7 + d' M+ a>D7,+ + d' M ant,- 



(88) 



(89) 
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with coefficients listed in Tables HVt IVHI 




71 > 



✓ 



K° 



FIG. 6: Feynman diagrams for K° — > 7r°7r° at one-loop ChPT. 



TABLE V: One loop contributions from dimension-7 operators (I) 



(L,R)ai 




K° — > Vacuum 




K 


b' K 




(15,3)i/ 2 ,s 





-12 


12 


(15,3)i/ 2lJ 4 


4 


-4 





(15,3)i /2 


-4 


4 





(6,3)i/2 


6 


12 


-18 


(3,15)i /2 ,l 


18 


-12 


-6 


(3, 15)l/ 2;J R 


6 


-12 


6 


(3, 6)i/ 2 ,L 


2 


4 


-6 


(3,6)i /2iR 


2 


4 


-6 



There are several comments we would like to make. First, the quark-mass-dependent 
operators (operators constructed with X±) contribute in the physical process K — > irir, in 
contrary to 2 in the SM case. The reason is similar to that for the higher-order operators 
of (27,1) (22)]: they cannot be expressed as a total divergence by the equations of motion. 
Similarly, these new operators will not act as a generator for rotation in s — d plane like the 

Therefore the one-loop matrix elements of K 
± irid) as they did at tree level. 



e 2 8,1) does 
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will no longer be 



31|, [33|, 

proportional to {m s 

Second, the masses appearing in our result are either bare masses or the renormalized one 
depending on the processes. For the unphysical processes K — > and K — > n, we use the 
bare masses, whereas for the physical K — ► tttt, the one- loop renormalized mass is implied. 
It makes the comparison to the experimental result feasible. 



Third, in 35 



the author claimed that infrared-sensitive terms like m 2 K log 
diverges in the m n —* limit, will emerge in the K — > irir matrix element. We have checked 



m^, which 
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TABLE VI: One loop contributions from dimension-7 operators (II) 

K° — > Vacuum 



(L, Ft) a t x+ 


K K 

77, A 


b' 

77, 7T 


b' K K 




b' K 


7r,7T 


(10, 0)l/2,S,X + 


A O / O 

—43/8 


27/8 


1 OO / A 

— 123/4 


111 I A 

111/4 


r 1 /o 

—51/8 


nn /o 

99/8 


(10, 0)l/2,S,X_ 


7/8 


L) 


—33/4 


U 


00 /o 

— 33/8 


—3 


(10, 6)i/2,A,X+ 


25/24 


— 17/24 


27/4 


1 cr / a 

— 15/4 


Ol 1 O 

—21/8 


/o 

—3/8 


(15, 6)l/2,A,X- 


— 25/24 


1 /Q 

1/3 


-3/4 


U 


—27/8 


U 


(15,3)i/2,x+ 


/o 

—0/0 


17/24 


1 r /,( 

— 15/4 


15/4 


/o 

—3/8 


/o 

3/8 


(15,3)i/2,x_ 


/o 

—3/8 


1 /Q 

-1/3 


1 K / A 

— 15/4 


U 


io 

—3/8 


U 


{ r o\ 

(6,3) 1/2jX+ 


/o 

3/2 


5/8 


—21/4 


9/4 


—57/8 


45/8 


(M)l/2,X- 


0/2 


1 /o 
1/2 


Ol 1 A 

— 21/4 


U 


—0 f /8 


9/ 2 


/"o 1 

(3, 15)i/2,L,X + 


25/8 


—5/8 


45/4 


OO 1 A 

— 33/4 


9/8 


a r /o 

—45/8 


(3, 15)i/ 2 ,L,X_ 


25/8 


1/2 


45/4 





9/8 


-3/2 


(3, lS)!^,/?,^ 


9/8 


-5/8 


33/4 


-33/4 


33/8 


-45/8 


(3, 15)i/2,k,x_ 


-9/8 


1/2 


-33/4 





-33/8 


-3/2 


(3,6)l/2,£,x4 


-43/24 


13/8 


-45/4 


33/4 


-9/8 


21/8 


(3,6)i/ 2 ,i,x_ 


-7/24 


1/2 


15/4 





3/8 


-3/2 


(3,6)i/ 2 ,K,x+ 


3/8 


-5/24 


3/4 


-3/4 


27/8 


-15/8 


(3, 6)i/ 2) H,x_ 


-3/8 


1/6 


-3/4 





-27/8 


3/2 



the result by keeping pion masses explicit in our calculations, and found that all such terms 
canceled when summing all the diagrams. Therefore it is safe to take the limit m n — > 0. 

Finally, there are a large number of unknown non-perturbative coefficients in the new 
operators. For dimension-6 operators, the traditional way of determining these coefficients 
by calculating simple processes like K + — > 7r + is suffice. In dimension-7 cases, however, the 
two simple processes K° — > and K + — > ir + are not enough to determine all the coefficients, 
unless there is the so-called CPS symmetry. Adding other simple processes like K° — > n and 
K° — > rj will not improve the situation since they are not independent in the SU(3) limit we 
are working on. We can in principle get more relationships when away from SU(3) chiral and 
isospin symmetries, but many more new coefficients will enter as well, and then we need even 
more relationships to determine all the coefficients. Therefore we could either rely on some 
model-dependent assumptions or calculate more complicate processes on lattice directly. In 
any case, the ChPT calculations can serve as a check for relations among coefficients from 
lattice or other non-perturbative model calculations. 

V. CONCLUSION 

The standard model calculations for the direct CP violation in non-leptonic kaon decay 
have not been entirely settled due to large cancelations between different matrix elements. 
Therefore, there is a considerable interest in understanding this phenomena from beyond 
standard model physics. However, we do not know yet what form the new physics will 
take, either supersymmetry, left-right symmetry, large extra dimensions, or little Higgs, or 
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TABLE VII: One loop contributions from dimension-7 operators (III) 



(L,R) A i 




K+ -» 7T+ 






K° -> 7r°7r° 




(4f) 




(c' M )x 


Ok) 




(d' K )x 


(15, 3)^/9 

V / u J of I 


208/3 


3/4 


-52/3 


8 


-2/9 


-2/9 


(15. 3 ) 1 /o g 


196/3 


28/3 


-49/3 


-640/9 


214/9 


160/9 


(15 3 li fn A 
V x ") Kj Jlf2,A 


-4/3 


8/3 


1/3 


-352/9 


70/9 


88/9 


(15, 3) 1/9 

V " / 1 f 2 


4/3 


8/3 


1/3 


-8 


2 


2 


(6. 3 )q /O 
V w > tJ Jof2 













8/9 _2/9 


-2/9 


(6,3)1/9 


-80 


12 


9 


-848/9 


212/9 


212/9 


(3, 15)3/2,1 


-16 


2/3 


26/3 











(3, 15) 3 /2,fl 


-16 


2/3 


26/3 


-16 


2/9 


-2/9 


(3,15)i/ 2 ,i 





-22/3 


29/3 





-3 


-3 


(3, 15) 1/2 ,_r. 


-32 


26/3 


23/3 


128 


-133/9 


133/9 


(3,6) 1/2 ,l 


-80/3 


-4 


3 


-8 


5 


-3 


(3, 6)i/ 2 ,B 


80/3 


4 


-3 





-97/9 


97/9 



something else. Presumably, the Large Hadron Collider will help us to identify it in the next 
few years. 

In this paper, we aim to study a general effective theory for non-leptonic kaon decay 
which has its origin from beyond SM physics. We systematically classify the dimension-5, 6 
and 7 quark and gluon operators according to their chiral structures. Using chiral symmetry, 
we derive tree-level relations between the matrix elements involving zero, one and two pions. 
This is useful because lattice calculations of multiparticle matrix elements are much harder 
than these for few particles. We have also calculated the leading chiral logarithmic behavior 
of these operators in ChPT. The result again will be useful for calculating matrix elements 
of these operators on lattice. We have not consider them in quenched QCD formulations, as 
the rapid progress in lattice QCD calculations makes queched studies much less useful than 
the past. 

We thank J. Bijnens and M. Golterman for useful correspondences. This work was 
partially supported by the U. S. Department of Energy via grant DE-FG02-93ER-40762. H. 
W. Ke acknowledges a scholarship support from China's Ministry of Education. 
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APPENDIX A: LEADING CHIRAL-LOGARITHMS IN SM OPERATORS 



The leading chiral- logarithms in SM operators have been calculated by many authors 
[H, El, 22, 23, 32], and for completeness we list the result here. Notice that the results 
quoted in Eqs. (80) and (81) in [36( contain sig n errors. The result for K° 



7r°7r° here is 

different from that in [l9|, as pointed out in |35j. Results for K° — > Vacuum and K + — > 7r + 
are presented in terms of bare masses and couplings, while for K° — > 7r°7r° we use the physical 
mass. 

The operators we use are defined in the main body of the paper (Eqs. [581- 165]) . We first 
consider the matrix elements between K° and the vacuum, 
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3 13 7 
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13 3 
- —L(m v ) - -L(m K ) - -L(m n ) 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 



where / is the bare meson decay constant, and m^^^rriKfi are bare masses of mesons. Due 
to the isospin conservation, only 1=1/2 part of the operator can contribute. 

For K + — > 7r + matrix elements, we apply a common mass tum for all the mesons. There- 
fore the momentum is conserved in the process. 
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(A6) 

(A7) 
(A8) 
(A9) 
(A10) 
(All) 



This result is useful in lattice calculations where the pion mass can be adjusted through 
quark mass parameters. The K° — > 7r° matrix elements can be obtained from the above by 
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using, 



(AOai=i/2\K°) = -^(7T + \0 AI=1/2 \K + ) 
(AOai=3/2\K°) = V2(tt + \O ai=3/2 \K + ) . 



(A12) 

Finally, for K — > -kit, we take the limit m n — > and keep the kaon mass dependency only, 
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1 - —L{m K ) 



(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 



(8 1) 

Here the physical mass of kaon is used. Note that the weak mass operator, 09 ' , will not 



contribute to the K — > nn matrix element as being pointed out in |30|, |31|, |33|, [34 . 

The matrix elements for the final state |7r + 7r~) is related to the above ones simply by 



A+- = -y=(A 2 + V2A ) . 
A 00 = J^(-V2A 2 + A C 



(A19) 



Compared with the angular momentum relation, the Aq amplitude has a factor of —a/2. 
The minus sign arises from the definition of w + = (n 1 + i-K 2 )/\/2 which has a different sign 
from the usual spherical tensor definition. The \/2 accounts for the identical particle nature 
of two 7r°'s, which is usually accounted from by a factor of 1/2 in the final state phase space. 
From the above relation, we derive: 



(7T + 7r-|O A/=1/2 |A ) = (7r 7r°|O A/=1/2 |A°) 



7T 7T 



\o 



A7=3/2 



K°) 



\(n°n \O AI=3/2 \K°) . 



(A20) 



Using the relation (1A12I) and (1A20I) . it is easy to check that the result for (8,8) operators is 



consistent with that in Ref. 23 
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